On the effective one-component description of highly asymmetric hard-sphere binary 

fluid mixtures 
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The phase diagram of a binary fluid mixture of highly asymmetric additive hard spheres is investi- 
gated. Demixing is analyzed from the exact low-density expansions of the thermodynamic properties 
of the mixture and compared with the fluid-fluid separation based on the effective one-component 
description. Differences in the results obtained from both approaches, which have been claimed to 
be equivalent, are pointed out and their possible origin is discussed. It is argued that to deal with 
these differences new theoretical approximations should be devised. 

PACS numbers: 05.70.Fh, 64.10.+h, 64.70.Ja, 64.75. +g 



It is well known that, due to the widely different time 
and length scales involved, it is a very difficult task to 
describe a soft matter complex fluid composed of meso- 
scopic (colloidal) particles and microscopic (small) par- 
ticles, i.e. a colloidal suspension. Therefore people have 
struggled for years trying to devise strategies that may 
render a tractable description. On the theoretical side, 
a rather successful and often employed approach in the 
colloids literature is that of coarse graining. In this tech- 
nique the idea is to integrate out the irrelevant degrees 
of freedom (the small particles) to end up with a one- 
component system of colloidal particles described by ef- 
fective interactions Q . One of the most celebrated effec- 
tive interactions results from the depletion mechanism 
first described by Asakura and Oosawa 0]- Suppose 
that two colloidal hard spheres are immersed in an ideal 
fluid of small hard spheres. When the separation of the 
surfaces of two large spheres is less than the diameter 
of the small ones, the depletion of the latter from the 
gap between the colloids leads to an unbalanced osmotic 
pressure which, in turn, results in an effective attraction 
between the two large spheres. Although the Asakura- 
Oosawa depletion pair potential is concerned only with 
the limit of extreme dilution, the theory of the depletion 
mechanism has recently been reexamined by also includ- 
ing the interactions between the small spheres HQ- In 
several of the recent studies the depletion mechanism has 
been invoked as the possible driving force for demixing, 
i.e. a phase transition whereby a binary fluid mixture of 
hard spheres separates into two fluid phases of different 
composition. 

The aim of this paper is to shed some more light on 
the usefulness or limitations of the depletion pair poten- 
tial picture, i.e. the consideration of an effective one- 
component fluid of large spheres with effective pair in- 
teractions for describing the phase behavior of a binary 
fluid mixture of highly asymmetric additive hard spheres 
H . Such an issue has recently been addressed by simula- 
tions H0H E| but the analysis unfortunately provides 
no definitive answer due to limitations in the presently 
available simulation algorithms. 



Recently jlOj we have analyzed the demixing transi- 
tion of a binary fluid mixture of N = Ni + N2 additive 
hard spheres of diameters a\ and a% (a± > (72)- For 
a fixed diameter ratio 7 = o^/fi, the thermodynamic 
properties of the mixture were described in terms of the 
number density p = N/V , with V the volume, the partial 
number fraction of the big spheres x = N\/N, and the 
total packing fraction 77 = 771 + 772, where 771 = irpafx/6 
and 772 = Trpal(l — x)j 3 /6 are the packing fractions of 
the large and the small spheres. Demixing was stud- 
ied by starting from the exact low-density expansions 
of the thermodynamic properties, using as input ana- 
lytical expressions II li Il2| and numerical computations 
111 13 III 13 IH Ell for the virial coefficients so far re- 
ported in the literature (up to the sixth). For all the size 
asymmetries considered we found that, by successively 
incorporating in the density expansions a new exact virial 
coefficient, the critical consolute point moves to higher 
values of the pressure and of the total packing fraction 
as shown in Fig. 1. It is also seen in this figure that the 
rate of convergence is slow and so definitive conclusions 
on the location of the critical consolute point cannot yet 
be reached with the available data on virial coefficients. 
We have also compared these findings with the fluid-fluid 
separation resulting from different empirical proposals for 
the equation of state 0, |3 , that attempt to repro- 
duce the virial behavior and/or to comply with consis- 
tency conditions of the contact values of the radial dis- 
tribution functions. We have realized that for historical 
reasons all these "rescaled" equations of state yield the 
exact second and third virial coefficients but inherit the 
singularity of the Percus-Yevick theory. We have shown 
that this apparently innocuous singularity arbitrarily in- 
troduces, as compared to the exact low-density expansion 
of the pressure, additional dependencies on x through the 
fourth, fifth, and so on virial coefficients. Such dependen- 
cies make the location of the critical consolute point to 
vary enormously from one equation of state to another, 
so that the predictions stemming from each of them are 
indeed unreliable. 

The demixing transition of a binary mixture of additive 
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FIG. 1: Critical total packing fraction r\ c as a function of the 
critical reduced pressure p* c , where p* = fipa\, (3 = 1/fcsT, 
with fcs Boltzmann's constant, and T the absolute temper- 
ature, as obtained from the exact low-density expansions for 
three values of the diameter ratio 7. Successive points (from 
left to right) indicate the addition of a new exact virial coef- 
ficient, starting with the second. Lines have been drawn to 
guide the eye. Dash-dotted line and triangles: 7 = 0.05; 
dashed line and squares: 7 = 0.1; solid line and circles: 
7 = 0.2. 



hard spheres may be analyzed using different thermody- 
namic planes. In order to compare to well documented 
results in the effective one-component fluid, we first con- 
sider the (771—772) plane. In Fig. 2 we present the binodals 
resulting from the exact low-density expansions up to the 
sixth virial coefficient for two size asymmetries, namely 
7 = 0.2 and 7 = 0.1, and up to the fifth virial coefficient 
for 7 = 0.05. Note that the density region spanned by 
each binodal cannot go beyond the point where the to- 
tal packing fraction occupied by the spheres reaches for 
the dense phase the highest possible value. Since these 
values are not known, we have drawn all the curves up 
to the limiting (unphysical) value 77 = 1. A comparison 
between these binodals and parallel ones derived using 
the depletion mechanism |f| indicates the following. For 
7 = 0.2, the low-density expansion predicts demixing, 
whereas there is no fluid-fluid separation in the effec- 
tive one-component fluid. The qualitative behavior of 
the binodals for 7 = 0.1 and 7 = 0.05 is, on the other 
hand, rather similar in both approaches. However, there 
are two main differences: 1) the values of 772 for the mix- 
ture are much lower than the corresponding ones in the 
effective one-component fluid; 2) the critical consolute 
points already lie outside the range where the binodals 
are predicted in the effective one-component fluid. 

Figure 3 is yet another presentation of the binodals, 
this time in the (771 — p) plane, with p the pressure. 
These binodals follow from the low-density expansions 
by keeping the same number of exact virial coefficients 
as in Fig. 2. Since no values of the pressure for the co- 
existing phases have been reported in the effective one- 
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FIG. 2: Binodals of asymmetric binary fluid mixtures of 
additive hard spheres in the (771-772) plane as obtained from the 
exact low-density expansions (continuous lines) for the same 
values of the diameter ratio as considered in Fig. 1. Filled 
circles indicate the critical consolute points. The dashed lines 
are the results of Ref. 0. The inset in the lower pannel is the 
curve for 7 = 0.05 in a magnified scale. 



component fluid in [f| , a direct comparison between both 
approaches is not possible in this case. Nevertheless, the 
following comments are in order. It is known that a one- 
component fluid with pairwise interactions consisting of 
a repulsive part and a short-range attraction can only 
produce in the (771 — p) plane a van der Waals-like fluid- 
fluid separation ending in an upper critical point, the 
concavity of the binodals pointing downwards. More- 
over, it is well established that, by strongly reducing the 
range of the attractions with respect to that of the re- 
pulsions, the usual phase diagram with two stable transi- 
tions, namely fluid-fluid and fluid-solid, changes and the 
fluid-fluid transition becomes metastable with respect to 
the fluid-solid transition [2^. By a further reduction of 
the range of the attractions, a stable isostructural solid- 
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FIG. 3: Binodals of asymmetric binary fluid mixtures of 
additive hard spheres in the (rji-p*) plane as obtained from 
the exact low-density expansions for the same values of the 
diameter ratio as considered in Figs. 0and|2] Dash-dotted 
line: 7 = 0.05; dashed line: 7 = 0.1; solid line: 7 = 0.2. Filled 
circles indicate the critical consolute points. 



solid transition develops at high densities in the phase 
diagram [2^, [24| . The phase behavior found in Ref . 
follows indeed this trend. Note that although the zero- 
body and one-body terms (the volume terms) in the ef- 
fective one-component description play no role in deter- 
mining the phase equilibria, they are however responsible 
for inverting the concavity of the binodals from down- 
wards to upwards, yielding the same concavity as the 
one found in the true binary mixture. The importance 
of the volume terms in the determination of the com- 
plete phase diagram should therefore not be ignored 0. 
On the other hand, the resulting phase diagrams found 
in Ref. [j| also indicate that, for the depletion pair po- 
tential used in this reference, the small repulsive barrier 
that follows the attractive potential close to the surface of 
the large sphere does not qualitatively modify the phase 
diagram that would result by only considering a purely 
short-range attractive potential. The predictions for the 
phase diagram in the effective one-component fluid |5j 
for 7 = 0.1 (a metastable fluid-fluid transition) and for 
7 = 0.05 (a metastable fluid-fluid transition and a stable 
solid-solid transition) are indeed at grips with the fore- 
going remarks. 

It should be clear that due to the fact that only a few 
exact virial coefficients are available and that the radius 
of convergence of the virial scries in the case of mixtures 
is not known, definite conclusions derived from the low- 
density expansions cannot be drawn. Nevertheless, the 
observed trends lead to the following considerations. As 
shown in Fig. 1, as one increases the number of terms 
in the low-density expansions, both the critical pressure 
and the critical total packing fraction grow. It is not clear 
whether this feature will change by including more exact 
virial coefficients (when available) in the density expan- 
sions, but all the evidence so far makes it unlikely. In any 



event, provided the trend obtained with the low-density 
expansions is a true feature, two possible scenarios arise. 
One may envisage the convergence of the critical con- 
stants to a thermodynamic state with either a physical 
or an unphysical critical total packing fraction. Indeed, 
the terms physical and unphysical depend of the resulting 
solid phase at high densities. For example, if the freez- 
ing transition corresponds to the partial freezing of the 
large spheres, i.e. a face-centered cubic lattice formed by 
the large spheres while the small spheres remain disor- 
dered [19| , the close packing fraction of the large spheres 
i]l P ~ 0.74 defines the boundary between physical and 
unphysical. In the first scenario the demixing transition 
would be thcrmodynamically metastable with respect to 
the freezing transition as found in previous studies, while 
in the second no demixing would be present at all. At this 
stage we have no way to decide which of the two scenarios 
is correct since, by incorporating up to the fifth (sixth) 
virial coefficient in the exact low-density expansions for 
7 = 0.05 (0.1), the critical packing fraction of the big 
spheres 77J is t}\ = 0.595(0.557). Concerning this point, a 
related analysis carried out with the Boublfk-Mansoori- 
Carnahan- Starling-Leland (BMCSL) equation of state 
[25] (which predicts no deximing) is illustrative. This is 
shown in Fig. The striking feature to be noted is the 
excellent agreement of the results obtained with the exact 
low-density expansion (c.f. Fig. and those obtained 
with the truncated BMCSL equation of state. By includ- 
ing additional virial coefficients of the BCMSL equation 
of state it is seen that r\ c ~ r\\ exceeds the close-packing 
fraction of the large spheres. 

In order to understand the differences between the 
phase diagrams presented in Fig. 2 using the exact low- 
density expansions and the effective one-component de- 
scription, the following considerations should be pointed 
out. For 7 = 0.2 both approaches would be equivalent 
whenever the critical packing fraction of the big spheres 
resulting from the low-density expansions converges to 
an unphysical value. As indicated above, for 7 = 0.05 
and 7 = 0.1 the critical consolute points obtained with 
a reduced number of exact virial coefficients already lie 
outside the range where the effective one-component de- 
scription locates the binodals. In these cases the effective 
one-component description would only provide a qualita- 
tive approximation of the true binary mixture if 77J con- 
verges to a physical value. Otherwise (convergence to 
an unphysical value) the predictions of both approaches 
would be completely at odds with each other in this 
basic issue. All these possibilities can indeed only be 
tested by including more exact virial coefficients in the 
density expansions. As a possible source for the differ- 
ences found between the exact density expansions and 
the effective one-component description we suggest that 
it may be due to the fact that, in order to determine in 
the latter the phase diagram and to perform the direct 
simulations |5j, the Carnahan-Starling equation of state 
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FIG. 4: Critical total packing fraction r/ c as a function of the 
critical reduced pressure p* as obtained from the truncation 
of the virial expansion of the BMCSL equation of state for 
two of the values of the diameter ratio 7 considered in Fig. Q 
Successive points (from left to right) indicate that an addi- 
tional (approximate) virial coefficient has been retained in the 
truncation, starting with the second. Lines have been drawn 
to guide the eye. Dash-dotted line and triangles: 7 = 0.05; 
dashed line and squares: 7 = 0.1. 



is used for describing both the big and the small hard 
spheres. Note that the Percus-Yevick-type singularity of 
the Carnahan-Starling equation of state may induce for 
each component different (through 771 and 772) spurious x- 
dependences similar to those discussed in Ref. [10| in con- 
nection with equations of state derived from "rescaled" 
virial expansions. These dependencies would in turn lead 
to unreliable predictions on the location of the demixing 
transition. 

Up to here and due to the involved approximations, 
none of the arguments and results that we have discussed 
lead to a neat conclusion concerning the claimed equiva- 
lence of the effective one-component description vis a vis 
the true binary mixture description of highly asymmetric 
hard spheres. Therefore, new theoretical approaches to 
deal with this problem should be sought. In this respect, 
systematic and controlled approximations in the fluid 
regime, may profitably be undertaken. Of course these 
are presently limited by the available virial coefficients. 
Also, the role of reliable simulations capable of exploring 
the region of interest can hardly be overemphasized. Ir- 
respective of the technical difficulties that no doubt are 
involved both in the simulations and in the computa- 
tion of additional virial coefficients of highly asymmetric 
binary fluid mixtures of additive hard spheres, we hope 
that the present paper may also serve to stimulate efforts 
in this direction. 
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